Rigorous derivation of single-atom responses from the triple scattering signal 
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We study coherent backscattering of intense laser light from three immobile two-level atoms using 
the master equation approach. The master equation is solved analytically, and the triple scattering 
spectrum is expressed in quadratures. We show that this solution can be obtained via a self- 
consistent combination of single-atom spectral responses, and is equivalent to the solution following 
from the diagrammatic pump-probe approach. We deduce the general expressions for single-atom 
building blocks which can be used in simulations of multiple inelastic scattering of laser light in the 
bulk atomic medium. 



I. INTRODUCTION 

Coherent backscattering (CBS) is an enhancement of 
the average intensity of the resonant wave backscattered 
off a dilute disordered medium due to the constructive 
interference between the counter-propagating multiply 
scattered amplitudes possessing time-reversal symmetry 
1 1 . Recent experiments on CBS of light from cold atoms 
2 -Q studied the impact of the atomic internal struc- 
ture on the phase coherence of the interfering waves. It 
was shown that the spin-flips of a photon on degenerate 
dipole transitions of rubidium atoms [HH Q , as well as 
nonlinear inelastic scattering from saturated strontium 
[j| and rubidium Q atoms induced by intense laser field 
significantly reduce the enhanced backscattering. 

Whereas there is a detailed description of CBS of light 
by atoms with arbitrary internal degeneracy in the elas- 
tic (linear) scattering regime 0, 0; @] > multiple inelastic 
scattering theory exists only for atoms weakly saturated 
by two photons [9] . But the parameter regimes of moder- 
ate or strong saturations probed in the experiments 0, H| 
requires a non-perturbative treatment of the atom-laser 
field interactions, and is beyond the reach of the diagram- 
matic approach put forward in Q . 

Standard non-perturbative quantum-optical methods, 
such as a master equation approach (see, for instance, 
[I3-[l2j]), can be used for a description of CBS of intense 
laser light from cold atoms. However, due to the expo- 
nential growth of the Hilbert space with increasing num- 
ber of scatterers these methods are fundamentally limited 
to systems of a few scatterers. So far, the master equation 
approach and extensions thereof have been successfully 
applied to treat double scattering [l3l - |l6l |. 

To avoid the above problem of the exponential growth 
of the Hilbert space, and yet to account for the atom- 
laser field interaction strength non-perturbatively, we 
proposed a stochastic method of solving the multiple 
scattering problem in dilute cold atomic clouds which 
we call the pump-probe approach to CBS [13, EH ■ The 
main idea of this method is to express the multiple scat- 
tering signal through single-atom responses to a classi- 
cal polychromatic random field, of which one component 
represents a laser (pump) field, while the remaining com- 
ponents describe weak (probe) fields scattered from the 
surrounding distant atoms. 



So far, the validity of the pump-probe approach was 
proven in the case of double scattering by two two-level 
atoms pj], EH by analytically establishing the equiva- 
lence with the results following from the master equa- 
tion approach [l9j]. Furthermore, we showed that the 
expressions for the double scattering background and 
interference spectra can be represented graphically, by 
self-consistently combining single-atom building blocks 
20]. Thus, the pump-probe approach is related to the 
non-perturbative methods of quantum optics, on the one 
hand, and to the diagrammatic scattering theories, on 
the other hand. 



In the present contribution, we explore these relations 
further. Namely, from the rigorous solutions of the mas- 
ter equation for three two-level atoms we deduce the an- 
alytical expressions for the triple scattering background 
and interference spectra of CBS. We present our results 
diagrammatically by self-consistently combining single- 
atom building blocks, analogously to the case of dou- 
ble scattering [20]. Thereby, we show the equivalence 
between the master equation and the pump-probe ap- 
proaches to CBS for triple scattering. Moreover, we de- 
duce general expressions describing the spectral response 
of an atom subjected to an arbitrary number of probe 
fields. These expressions will be required in future work 
to implement the pump-probe approach for a medium 
consisting of an arbitrarily large number of atoms. 



The structure of this paper is as follows. In the next 
section, we present the model, the master equation, and 
the main quantity of interest - the triple scattering con- 
tribution to CBS of light from three two-level atoms. In 
Sec. IIIII we describe the method to find an analytical so- 
lution of the master equation which links it to the pump- 
probe approach. Sections IIVI and |V] contain our main 
results. They include diagrammatic expansions of single- 
atom building blocks into the elastic and inelastic spec- 
tral responses whose composition yields the ladder and 
crossed spectra (Sec. IIV[) . and the general expressions for 
single-atom building blocks (Sec.|V|. We present analyt- 
ical and numerical triple scattering spectra in Sec. I VII 
Conclusions of this work are presented in Sec. I VIII 
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FIG. 1: Three atoms (gray circles) located in the radiation 
zone of each other are driven by a quasi-resonant laser field 
of arbitrary intensity with a wave- vector k^. They radia- 
tively decay (dotted arrows) into a common electromagnetic 
bath, which also mediates the resonant dipole-dipole inter- 
action (double wavy arrows). We will be interested in the 
far-field intensity of the scattered light along k. 



II. MODEL 



inversely proportional to the distance between the atoms 
A and /i: 

7V^- 7 (l-(d.f v ) 2 )- , (2) 

where d and f \^ are unit vectors along the atomic dipole 
moment and radius-vector connecting the atoms A and 
/i, respectively, and = \v\ — r^|. It should be stressed 
that Eq. ^ is valid only in the limit k^rx^ 1, which 
is standard in the context of CBS, and is also assumed 
here. Note that although the constants T Xlx = T^ x and 
T x = T* x , the interaction Liouvillians in Eq. ([TJ asso- 
ciated with these constants are not symmetric under the 
permutation of the indices A and /i. Physically, they de- 
scribe different excitation transfer processes induced by 
the far-field dipole-dipole interaction Liouvillian . We 
will discuss this issue in Sec. MI C II 



A. Master equation 

We will start out with introducing our model depicted 
in Fig. [TJ It includes three quantum scatterers (two-level 
atoms) embedded in a common electromagnetic reser- 
voir at fixed random positions r\ and excited by a quasi- 
resonant continuous wave laser field with a wave vector 
k^. We will be interested in the disorder-averaged far- 
field intensity of the light scattered into a given direction 
k. 

To describe CBS of light from three laser-driven atoms, 
we will employ the master equation approach [ToL [l2| 
which has already been used for calculating double scat- 
tering from two saturated Sr [l3|, [l6| and two- level atoms 
19]. Therefore, we will next present a master equation 
governing the evolution of the quantum mechanical ex- 
pectation values of an arbitrary atomic observable Q of 
a three-atom system. In the Heisenberg picture and in 
the frame rotating at the laser frequency u>l, (Q) obeys 
the following equation of motion [19| : 

3 
A=l 

- 1 (a+a x Q + Qa+a x -2a+Qa x )) (1) 

3 

A^=l 

Here, a x = |1) A (2| A and a+ = |2) A (1| A , with |1) A and 
1 2) a being respectively the ground and excited states of 
atom A, denote the atomic lowering and raising opera- 
tors. Furthermore, il\ = f2e lkl/ ' rA is the Rabi frequency 
dependent on the atomic position Y\, S = lol — loq is the 
laser-atom detuning, and 7 is half the radiative decay rate 
of the excited state. The lower line of Eq. ([1]) describes 
the retarded dipole-dipole interaction proportional to the 
complex couplings T\^. The couplings 7a m , in turn, are 



B. Triple scattering intensity 

Solution of Eq. ([TJ gives access to the expectation val- 
ues of observables of the field radiated by a system of 
three atoms. We will focus on the triple scattering con- 
tribution to the stationary scattered light intensity in the 
direction of the wave vector k. This quantity is a part of 
the general expression for the intensity of the scattered 
light given by 

(j(k)> M =5>A*r>-+ j2 k^ss^-^h (3) 

A A,^A 

where the angular brackets (. . .) ss denote the quantum- 
mechanical average in the steady state. The first and 
second terms in Eq. ([3]) describe the incoherent (or back- 
ground) and the interference intensities, respectively. 

In the framework of the master equation approach, the 
multiple scattering contributions correspond to the sub- 
sequent terms in a perturbative expansion of the sta- 
tionary solution of Eq. (UJ in the power series of the 
dipole-dipole interaction constants. The first-order term, 
proportional to \T Xfl \, describes the amplitude of double 
scattering between the atoms A and /x. Consequently, 
the intensity of triple scattering between atoms 1, 2 and 
3 must be proportional to IT12T23I 2 . Furthermore, the 
CBS signal is observed after the configuration, or disor- 
der, averaging over the random positions of the scatter- 
ers, which will be denoted as (. . .) CO nf ■ Summarizing the 
above, we can write the expression for the triple scatter- 
ing contribution to CBS as follows: 

«J(k))W)conf = E^A^Uf (4) 
A=l 

+ £ ((4®^)We ik -^- r -)) conf , 

A,^A 
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where the superscript (4) denotes the fourth order in the 
dipole-dipole coupling. Following the standard nomen- 
clature in the field of coherent quantum transport (2lj . 
we will refer to the configuration averaged background 
and interference contributions as the ladder and crossed 
intensities, respectively. 

In the following, we will ignore recurrent scattering, 
which is proportional to |I\ M | 4 . Although in a three- atom 
master equation, the recurrent and triple scattering con- 
tributions are of the same order of magnitude, in a dilute 
atomic cloud, the impact of the recurrent scattering on 
the observed signal is by far dominated by that of the 
multiple scattering. 

We will identify all triple scattering contributions 
which survive the disorder averaging. Since some of them 
are equal to each other, it will be sufficient to find only 
the non-equivalent contributions. To this end, in the next 
section we will study the structure and solutions of the 
master equation ([T]) in more detail. 



The steady state solution of Eq. © which is of the 
fourth order in the matrix V reads: 



(Q) 



(4) 



(GV^GA, 



(10) 



where G = —A l . Given the matrices A, V, and A, it 

is easy to find (Q)it^ numerically for any random posi- 
tions Yi of the atoms. Averaging this result over disorder 
gives ((Q^ 4 -')ss)conf and, consequently, the terms from the 
right hand side of Eq. (|4| contributing to the ladder and 
crossed intensities. However, in this work we would like 
to prove the equivalence of the results of the master equa- 
tion and the pump-probe approaches. Therefore, we will 
adhere to the analytical tools developed in [ill [2(3] . 

We will proceed with discussing the structure of the 
evolution matrices and presenting the recurrence rela- 
tions connecting the steady-state vectors x, y, and z. 



B. Recurrence relations 



III. SOLUTION OF THE MASTER EQUATION 

A. Basis set of operators 

Let us consider the evolution of the expectation value 
of an operator Q from the complete three-atom basis set 
of operators: Q G {<fi <g) q 2 <8> 9*3}, with 

qx = (l x ,<T-,a+,<T z x ) T , 

!a = cr+ a' + a{ = a+a^ - a^a+. (5) 

The normalization condition Tr [li <S> I2 <H> I3P] = 1 which 
holds for any atomic density operator p reduces the size 
of the operator space by one, leading to a set of 63 opera- 
tors. Their expectation values will be ordered as elements 
of the vector 



(Q) = (x,y,z) q 



(6) 



with 



x=((a 3 ),(a 2 ),(a 1 )) T , (7a) 
y = ((a 2 ® a 3 ), (ai (g) a 3 ) , (a 1 ® a 2 )) T , (7b) 
z = (a 1 <g> a 2 ® a 3 ), (7c) 

where (a\) is a Bloch vector corresponding to atom A: 

(^a) = (K-),K),K)) T . (8) 

The evolution of the vector (Q) is governed by the 
equation of motion: 



(Q) = (A + V)(Q) + A, 



(9) 



where the matrices A and V govern the dynamics of in- 
dependent and dipole-dipole interacting atoms, respec- 
tively, and the vector A is given by Eq. (|A4I) in Appendix 

m 



For a vector (Q) given by Eq. ©, the matrices A and 
V have the following block structure (see Appendix [A"|l : 

M+ \ / Un 

.1=1 L + My , V = \ U L U x VFn I . (11) 
L x N x J \ W x 

A perturbative expansion of the steady-state solution of 
Eq. ^ in the power series of the matrix V yields a system 
of coupled recurrence relations for the vectors x^ n > , y ( n > , 
and z 

x W =G+C/ny( n - 1 \ (12a) 

=G x U^x ^ +G x U x y^ +G x L+x W 

+ G x W^z( n ~ 1 \ (12b) 

z (") = D x W L y {n - 1] +D x W x z ( - n ~ 1) +D x L x yW 

(12c) 

with the initial conditions determined by (Q)is = for 
n < 0, and (Q)$ = (x (0 \ y (0 \ z (°)) T , where 



' (0) = ((o ; 3 ) (0) ,( ( ? 2 ) (0) ,(a 1 )^) T , 



(13a) 



y (0) = ((a 2 ) (0) ®( O ; 3 ) (0) ,( O ; i) (0) ®( ( ?3) (0) , 

( O ; i) (0) ®( O ; 2 ) (0) ) T , (13b) 
= (a 1 )W®(a 2 )W®(a 3 )(°\ (13c) 

and (cta) (0) = G\L. In Eqs. ([12]) and (JTSJ), we introduced 
a shorthand notation G\ = G\(0) (A = 1,2,3), G a = 
G a (0) (a = +, x), D x = D x (0), where 



G x (z) 



z-Mx 



G a (z) 



z-M, 



' , D x (z) ' 



z-N x > 
(14) 



are the Green's matrices. 
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V, 



> 4 . 



(a) 



(b) 



(O 



(d) 



Atom X 



V 



Atom n 



FIG. 2: Elementary excitation transfer processes between 
atoms A and fi generated by the dipole-dipole interaction ma- 
trix V. Black spots correspond to the atomic states after the 
exchange process has occurred, and solid (dashed) arrows - to 
positive (negative) frequency amplitudes. The positive (neg- 
ative) frequency amplitude of a transfer process from atom A 
to atom /j, is described by the matrix V^x (d) (Vx^ (a)), while 
the reversed process is described by matrix Vj^ (b) (V^a (c)). 



The possibility of demonstrating the equivalence be- 
tween the results of the master equation and the pump- 
probe approach is based on the analytical solution of the 
recurrence relations (|12p for fixed atomic positions, and 
on the subsequent analytical configuration averaging pro- 
cedure. 

In the following, we will identify the triple scattering 
paths for which configuration averaging of x ^ and y ^ 
yields the non-vanishing contributions around the back- 
wards direction k = — k^. 



C. Selection of triple scattering paths 

1. Interaction matrices and amplitudes of excitation tranfer 

In the framework of the master equation approach, 
multiple scattering paths can be defined according to the 
physical meaning of the interaction matrix V as describ- 
ing the excitation transfer processes between the atoms 
[1 91 ] . This matrix can be decomposed into a sum 



V 



3 

E 



(15) 



with V\n oc T\p and oc T£ . Each term in the right 
hand side of Eq. (|15[) has the same block structure as 
the matrix V itself (that is, given by Eq. (fTTj) ). and gen- 
erates an amplitude for a particular excitation transfer 
process between the atoms A and \i (see Fig. [2]). This 
amplitude corresponds to a propagation of the positive- 



(solid arrow) or ne gati ve- frequency (dashed arrow) field 
between the atoms [19| . Accordingly, the block matrices 
(C/jJ ) a , (W%J Q (a = n , l, x) are associated with a prop- 
agation of the positive-, and {U\^) a , (W\^) a with the 
negative-frequency amplitude between the atoms A and 

After the disorder averaging, those triple scattering 
processes which can be reduced to each other by relabel- 
ing the atomic indices result in identical contributions. 
Therefore, it is instructive to identify distinct types of 
triple scattering paths and consider one representative 
thereof. The intensity (0| then follows after multiply- 
ing the expression corresponding to a given type of triple 
scattering process by the number of times it occurs. 

Without restricting the generality, subsequently we 
will be concerned with the triple scattering processes 
whose positive frequency amplitudes include the path 
ri — > r 2 , and are thus proportional to T 2 \. It 1S easy to 
see that there are two triple scattering positive-frequency 
amplitudes containing the scattering process ri — > r 2 : 
**i — > T2 — )• r3 and ri — > r 2 r3. The first and sec- 
ond amplitudes are proportional to T 21 T^ 2 and T 21 T 23 , 
respectively. 

Because the atomic positions are random, the dipole- 
dipole coupling constants (see Eq. ^) carry random 
phases e tkLrXfi , for which we assume, since k^rx^ ^> 1 
as stated above, a uniform distribution between and 
2tt. Then, to survive the disorder averaging, the positive- 
frequency triple scattering amplitudes must be combined 
with corresponding negative-frequency (conjugate) am- 
plitudes. On top of that, the interference terms in Eq. (|4]) 
have additional, k-sensitive random phases. This im- 
poses further restrictions on the conjugate amplitudes. 

For a given set of four interaction matrices describing 
a particular triple scattering process there are 4! = 24 
different terms (due to the non-commutativity of the ma- 
trices). Fortunately, the phase of each of these 24 terms 
is the same. Therefore, to see whether a given process 
survives the configuration averaging it suffices to consider 
only one term. 

Using the above restrictions and analyzing different 
triple scattering processes, we have identified four non- 
equivalent contributions which survive the configuration 
averaging: two ladder and two crossed contributions. Be- 
low, we will discuss them separately. 



2. Ladder intensity 

The ladder intensity is given by the first line in the 
right hand side of Eq. We established that either 
of the two positive-frequency triple scattering processes, 
ri — s- r 2 — > r 3 and ri — > r 2 r3, can be combined with 
one negative-frequency triple scattering process, to yield 
a phase-independent contribution. 

For the first excitation transfer process, proportional to 
T 21 T 3 * 2 and generated by matrices V 21 , V£ 2 , the conjugate 
amplitude is generated by matrices Vi 2 , V23, and leads 
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to a phase- independent solution for {a^a^)H\ The cor- 
responding triple scattering process which is composed 
of the co-propagating complex conjugate amplitudes, is 
shown in Fig. [Ha). Here, atom 3 is the source of the 
backscattered signal, depicted by one dashed and one 
solid arrows which are associated with the operators a + 
and o~ , respectively. 

By simply exchanging the atomic indices 1, 2, 3, it 
is easy to see that the total number of triple scattering 
processes which yield the equivalent expressions upon the 
disorder averaging is equal to 3! = 6. 

Analyzing the second type of the triple scattering pro- 
cess which is generated by matrices V 2 *i, V 2 %, we iden- 
tified that, the only conjugate process which leads to a 
contribution to the ladder intensity is generated by ma- 
trices V\2 and V32. This triple scattering process results 
in a phase- independent solution for (cr 2 a 2 )^ > an d is de- 
picted in Fig. [3(b). Since the corresponding term is sym- 
metric with respect to a permutation between the indices 
of the outside atoms, there are in total 3 contributions of 
this type. 

Combining the contributions of both types, which we 
call "ladder type 1" and "ladder type 2" in the following, 
we obtain the following expression for the total ladder 
intensity 



(3) 



(3) 



(3) 



k, 



k % — 

(a) 





k, 



k 

(c) 



O r 3 



(b) 



4 - ' ' *A 



k 

(d) 



O h 



FIG. 3: Triple scattering processes surviving disorder aver- 
aging. Co-propagating amplitudes in diagrams (a) and (b), 
called "ladder type 1" and "ladder type 2", in the following, 
have no overall phase and contribute to the background inten- 
sity; (c) interference of counter-propagating amplitudes be- 
tween atoms 1 and 3 ("crossed type 1"), and (d) interference 
of amplitudes resulting from reversal of the amplitudes con- 
necting atoms 2 and 3, carry phase factors exp{j(k + k£) ^13} 
and exp{i(k + ki) • r23} ("crossed type 2"), respectively (see 
Eq. (1211) and preceding discussion), and thus contribute only 
in backscattering direction k = — k_L where these phases van- 
ish. 



(16) 



with 



3. Crossed intensity 



(17) 

where, from now on, we will for brevity drop the subscript 
'ss' in quantum mechanical steady state averages. Fur- 
thermore, the total ladder intensity splits into an elastic 
and an inelastic component: 



(3) 
J tot,l 



-(3) 



-(3) 



-(3) 



(3) 

1,1 T ^incl.l' -^tot.2 — ^el.2 ^incl.2' 



(3) 



(18) 



where the elastic component results from factorizing the 
expectation value in Eq. JT7J, ke., 



i2 ) i=6«o(a 3 -)W) COIlf 



(4)\ 



(19) 



(3) 

and similarly for Ly,. The non-factorizable inelastic re- 
mainder i = ^tot.i - ^el i> * — 1; 2, finally, is obtained 
as an integral over the frequency v of the detected pho- 
ton: 



(3) 



(20) 



(3) 

where the spectrum L\J Q \i{v) results from the Fourier 
transform of the non-factorizable part of the correspond- 
ing atomic correlation function (er + (i+T)er~ (t)) — (<r + (t+ 
T))(cr~(t)) occuring in Eq. (jTT)) with respect to r. 



Let us start discussing the crossed intensity by consid- 
ering the processes oc T 21 T% 2 . In accordance with our 
interpretation of the excitation transfer processes and 
their association with the atomic raising and lowering 
operators, atom 3 now emits a positive frequency ampli- 
tude contributing to the CBS signal (solid arrow). Cor- 
respondingly, the negative frequency amplitudes (dashed 
arrows) of the outgoing field can be emitted either by 
atoms 1 or 2. Therefore, there are two types of two-atom 
correlation functions, (af ® cif) an d {<J% ® a a) 1 
which contribute to the interference intensity, and will 
be called "crossed type 1" and "crossed type 2" in the 
following. 

As follows from Eq. the atomic dipole correlation 
function (a^ (g> cr^}^ (type 1) contributes to the CBS 
signal upon the configuration averaging if it includes the 
position-dependent phase factor which cancels itself with 
the phase exp(ik • r^). We found that this happens only 
for the counter-propagating amplitudes which are pro- 
portional to T 2 * 1 T 3 * 2 T2iT32 (Fig. 0(c)). In this case, the 
function (erf ® ^3)^ carries the phase exp(ik^ • 1*13) 
and, hence, the interference term survives the configura- 
tion averaging in the exact backscattering direction (that 
is, for k = — ki [22[)- A contribution of such counter- 
propagating amplitudes to the backscattered signal is 
generic for the CBS effect [H|. 

The term {a 2 <8> c r 3 ~)( 4 - ) (type 2) leads to another type 
of crossed contribution. It survives the disorder averag- 
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ing in the backwards direction for the scattering process 
oc T21T32T12T32. For this contribution, the negative fre- 
quency amplitude connecting atoms 2 and 3 is reversed 
(see Fig.[3fd)), whereas the one connecting atoms 1 and 2 
is co-propagating with the positive-frequency amplitude. 
Such a type of interference processes is characteristic also 
for CBS from nonlinear classical scatterers foL [23l|. 

Last but no least, there is an interference contribu- 
tion whose positive-frequency amplitude is proportional 
to TJjT^. The diagram for the corresponding triple scat- 
tering process which survives the disorder average can be 
obtained from the diagram in Fig. [3Jd) by substituting 
each solid arrow by a dashed one and vice versa. As the 
corresponding correlation function, (a^ ® ) " , is the 
complex conjugate of the expression corresponding to a 
diagram in Fig. Eld), its contribution can be accounted 
for by taking twice the real part of the expression de- 
scribing Fig. [Hd). 

Finally, by noting that the degeneracy of each of the 
crossed contributions is equal to 3! = 6, we can write 



r (3) _ r (3) 
°tot — '"'tot,! 



a 



(3) 

tot, 2' 



with 



C t ( l 2 =12Re[<<a+®a 3 -)( 4 V k - r -) 



conf j 



(21) 

(22) 
(23) 



Like in the ladder case, also the crossed intensity splits 
into an elastic and inelastic component C f 



(3) 



.(3) 



tot.i el,i 



(3) 

Cj no |-, i — 1,2, defined as the factorizable and non- 
factorizable components of the atomic correlation func- 
tion. 



IV. DIAGRAMMATIC REPRESENTATION OF 
THE TRIPLE SCATTERING CONTRIBUTION 

Our goal is to obtain analytical expressions for the 
triple scattering ladder and crossed intensities defined 
in the previous Sec. IIII C[ and to present them as com- 
pounds of single-atom spectral response functions. Since 
the derivation of these single-atom response functions 
starting from the master equation amounts to a quite 
lengthy calculation, we will, for the sake of clarity, first 
introduce the diagrammatic representation of the single- 
atom response functions and define the rules accord- 
ing to which these single-atom building blocks are con- 
nected with each other, before showing that the thereby 
diagrammatically obtained triple scattering signals are 
equivalent to the expressions derived from the master 
equation. 



A. Triple vs. double scattering 

For this purpose, let us recall the case of double scat- 
tering. Here, a diagrammatic representation in terms of 



Double scattering 

(A) 

O 





(b) 



FIG. 4: Single-atom blocks contributing to double scattering, 
(a) Background and (b) interference contributions. Dashed 
frames subdivide multiple scattering contributions into single- 
atom building blocks. 



sing le-atom building blocks has already been introduced 
[ll GJ, [20[ and been proven to be equivalent with the 
double scattering signal derived from the master equa- 
tion [ll|. We will therefore use this insight in order to 
construct analogous diagrams for the case of triple scat- 
tering. 



1. Multiple scattering as a combination of single- atom 
responses 

To start with, we plot the fundamental processes which 
survive the disorder average for the double and triple 
scattering contributions in Figs. |4] and El respectively. 

Let us consider the diagrams in Fig. |4] which depict 
the configuration averaged (a) background and (b) in- 
terference contributions for double scattering. Here and 
henceforth the laser field is not displayed, in order to 
lighten the diagrams. As shown in [17|, [lg], both the 
ladder and crossed spectra can be obtained by combin- 
ing single-atom building blocks which are enclosed in 
dashed frames (A)-(D) in Fig. B|a), (b). The interac- 
tion between the atoms is accounted for perturbatively 
via the fictitious classical probe fields (solid and dashed 
arrows) which are connecting the atoms. These ficti- 
tious fields represent the far-field dipole-dipole interac- 
tion between the atoms ff^|. Building block (A) is of 
the zeroth order, blocks (C), (D) of first, and block (B) 
of second order in the probe field or, equivalently, the 
dipole-dipole coupling strength. Figure [5] shows the com- 
position of the configuration averaged triple scattering 
diagrams in Fig. [5] using single-atom building blocks en- 
closed in dashed frames. Comparing the dashed frames 
in Figs. 0] and [S] we impose the compatibility condi- 
tions: (A')=(A")=(A), (B')=(B")=(B), (C')=(C), and 
(D')=(D). 

This simple analysis allows us to fully specify the struc- 
ture of the contribution to the ladder spectrum described 
by Fig. [SJa). It can be compounded from the building 
blocks (A) and (B) which can be found by solving the 
optical Bloch equations (OBE) under a classical bichro- 
matic driving [13, [H| . 

The second ladder contribution, Fig. \5[b), as well as 
the crossed diagrams, Fig.[5jc) and (d), contain three new 
building blocks, (E), (F), and (G), respectively. These 
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Triple scattering 





FIG. 5: Single-atom blocks contributing to triple scattering, 
(a) and (b) background; (c) and (d) interference contribu- 
tions (corresponding to Fig. EJa), (b), (c), and (d), respec- 
tively). Dashed frames subdivide multiple scattering contri- 
butions into single-atom building blocks. 



(a) 

= <or)W 

(b) 
(0 

= <a-(uu)> (+) 
(d) 

^O— » = <o-(u),oo)} ( - +) 



V-UJ 

(-) <■ - - 



(e) 

JJ 

^M-^> ~ P H (w;v) 

(g) 

v 

±||---> = P< + >(w;v) 

(h) 



- - -> w& ■ - -> 

UJ v 



= P <+ '(uj.cjoiv) 



FIG. 6: Full set of elementary single-atom building blocks 
for the double scattering contribution. Circles are (hatched) 
blank if the frequency of the outgoing arrow is (different 
from) equal to the laser frequency, (a) Complex scatter- 
ing amplitudes expressed through the induced atomic dipole 
moment and (b)-(d) corrections thereof due to the incom- 
ing probe fields; (e)-(f) Building blocks associated with the 
atomic dipole correlation functions describing inelastic scat- 
tering. In the response functions on the right hand side, the 
superscripts (0), (+), ( — ) and (+ — ) or ( — h) correspond to 
the number and type (solid or dashed) of incoming arrows. 



blocks describe the spectral response of the intermediate 
atom, which receives probe fields from its two neighbors. 
Therefore, within the pump-probe approach the blocks 
(E), (F), and (G) can be determined by solving the OBE 
under a classical trichromatic driving field. 

Because the blocks (A), (B), (C), and (D) familiar from 
the treatment of double scattering are important also for 
triple scattering, we will next recall the corresponding 
single-atom spectral responses. 

2. Diagrammatic pump-probe approach for double 
scattering 

It is convenient to represent the solutions of the OBE 
under a bichromatic driving from which the single-atom 
building blocks (A), (B), (C_h and (D) (see Fig. gj are 
deduced diagrammatically [20], as displayed in Fig. [5] 

Any elastic and inelastic spectral response function de- 
scribing the building blocks (A)-(D) can be compounded 
from circles and boxes shown in Fig. [5] Circles are asso- 
ciated with expectation values of the atomic dipole mo- 
ment, and boxes with Laplace transforms of the steady- 
state atomic dipole temporal correlation functions [20j. 
Both types of the averages are evaluated perturbatively 
to second order in the weak probe field, which is repre- 
sented by the incoming arrows in Fig. |51 The correspond- 
ing mathematical expressions are given in Sec. |V] 

Hatching the boxes and some of the circles indicates 
that the outgoing arrows have non-zero detunings from 
the laser frequency. This helps us to identify the fre- 
quencies of the outgoing arrows in the composition of 



the double scattering diagrams shown in Fig. HJa),(b) - 
by virtue of the physical interpretation of the elementary 
building blocks as the effective nonlinear susceptibilities 
[2(i| which was partially inspired by the earlier work of 

Moiiow dHH. 

In all diagrams in Fig. [6j the important information 
about the incoming and outgoing arrows is encoded in 
their type (dashed or solid), direction (incoming or out- 
going) , and frequency (more precisely, the detuning from 
the laser frequency). By convention, we do not label 
the arrows that are on-resonant with the laser frequency. 
The notation of the spectral response functions on the 
right hand side of each diagram in Fig. [5] embodies the 
relevant information about the shapes and arrows of the 
elementary blocks, and is explained in detail in Sec. |Vj 

In general, the diagrams in Fig. [5] correspond to the 
complex- valued spectral response functions. Their conju- 
gate diagrams are obtained by replacements of each solid 
arrow by a dashed one and vice versa. To obtain the cor- 
responding spectral response functions, we must flip all 
signs in the symbolic definitions thereof. For instance, 
«a-( W ,c))(-+))* = (a+(w,w)>(+-> = <(7+(w,w)>C-+>, 
where the latter equality follows because the incoming 
arrows are frequency-degenerate. By the same token, 
(P(+)(w;i/))* = P(-)(u;i/), whereas P (+ ~\tjj,uj;v) = 
P( ^'(u),u); v) is a real function. 

How to compose the building blocks (A)-(D) from the 
elementary diagrams shown in Fig. [6] was discussed in 
Ref. [2(j. The same setting will be generalized below 
in Sec. IIVI to describe spectral response functions with 
arbitrary number of the probe fields. 
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B. Ladder contribution 



Type 1 



(A) 

OTTt = 



(al) 

O— ► 
x 

o— > 



u2) 



The terms yielding the ladder spectrum of type 1 are 
shown in Fig. [7] in the form of three graphical equations. 
The left hand sides of these equations depict the single- 
atom building blocks (A), (B), (B) from Fig. BJa). To 
allow for an easy identification of the fields contributing 
to the detected signal, here and henceforth the arrows 
representing the backscattered fields are directed down- 
wards. Apart from this distinction, which becomes im- 
portant when combining the single-atom responses into 
triple scattering contributions, the two building blocks 
(B) in Fig. [7] are equivalent. 

The right hand sides of the graphical equations from 
Fig.[7]are the expansions of the building blocks into sums 
of the elastic and inelastic spectral response functions 
shown in Fig. [51 The arrows are not labeled by the fre- 
quency values in Fig. because these values are uniquely 
defined at the stage of self-consistent combination of the 
single-atom responses into triple scattering diagrams (see 
below Sec. UVDl) . 

Each of the elastic responses is depicted as a product 
(denoted by ' x ') of two elementary building blocks from 
Fig. |6ja-d) (as well as their complex conjugates). There 
are 2 n such products - the number of ways in which n in- 
coming arrows can be distributed among the two circles. 
Last but not least, every building block contains one el- 
ementary inelastic building block depicted by a hatched 
box. The boxes (a2) and (b5) in Fig. [7] correspond to 
diagrams of Fig. M, e ) and (h), respectively . 

Thus, the block (A) is equal to a sum of 2 terms, 
whereas the block (B) is equal to a sum of 5 terms. Con- 
sequently, a compound (A)(B)(B) yields the expression 

for the ladder intensity i of type 1, see Eq. ([TT]). con- 
sisting of 50 terms. Note that, as discussed in Sec. lIII C~2l 
see Eq. ([T7[) . the sum of these terms must be multiplied 
with a factor 6 in order to account for all different possi- 
bilities of exchanging the three atoms with each other. 



2. Type 2 



(B) 



(bi) 

o---> + 

x 



(b4) 



(b2) 

to-> 



(b5) 



(b3) 
X 

>o > 



(B) 

to 



-*o---> + 

x 

->o^ 



(bl) 



(b2) 



tO*0 + TT^OxO + 



•o 



4> v 

(b4) 



; i --->o x ; 



>o + 



->o x * 



(b5) 



FIG. 7: (Left) single-atom building blocks contributing to 
the ladder spectrum of type 1; (right) expansion of the 
single-atom building blocks in elementary single-atom build- 
ing blocks (see Fig. [6]). Circles and boxes describe the elastic 
and inelastic responses, respectively. The backscattered sig- 
nal is depicted by the downward arrows. 



originate from different atoms. Therefore, in general, the 
values of the incoming frequencies lo\ and U2 are not 
equal to each other. In the framework of the pump-probe 
approach, the evaluation of the corresponding spectral 
response functions requires solving single-atom OBE un- 
der a classical trichromatic driving. In the degenerate 
case wi = LJ2 = the functions {a^ (011,0)2)) ( h ' are re- 
duced to the functions (cr ± (w , w))*- ^ familiar from dou- 
ble scattering (see Fig. [6j). 



To obtain this contribution, we perform a diagram- 
matic expansion of Fig. [5jb) into elementary single-atom 
building blocks, as shown in Fig. [5] Since the middle 
atom receives four arrows, the decomposition of the elas- 
tic response functions consists of 2 4 = 16 terms, see 
Fig.[i;el-el6). 

Obviously, some of the elementary blocks in Fig. [Side- 
scribing the response of the middle atom, in particular, 
those with more than 2 incoming arrows and some with 
2 incoming arrows, are not present in Fig. [6J Therefore, 
we depict in Fig. [S] the new spectral response functions 
required for describing the building block (E). All ele- 
mentary blocks in Fig. [§] contain incoming arrows which 



Now, all elementary building blocks in Fig.[5]can be ex- 
pressed in terms of the elementary building blocks shown 
in Figs. [5] and [3] or in terms of their complex conjugate 
diagrams. Again, the mathematical expressions for the 
elementary building blocks shown in Fig. [H] (as well as 
those shown in Fig.[TTJ see Sec. UVCj) are given in SeclVl 
see Eqs. (H7J and (joT)) . 

By combining all diagrams (A), (E), (A), we obtain 
68 = 2x17x2 terms, the sum of which must finally be 
multiplied by a factor 3 in order to obtain the expression 
-ktot 2 for the ladder intensity of type 2, see Eq. (TlT[) . 
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FIG. 9: Elementary building blocks required to describe the 
response of the intermediate atom in Fig. [8] account for 2, 3, 
and 4 arrows with frequencies uii and u)2- (a)-(g) scattering 
amplitudes; (h) frequency correlation function. 
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FIG. 8: (Left) single-atom building blocks contributing to 
the ladder spectrum of type 2; (right) expansion of the 
single-atom building blocks in elementary single-atom build- 
ing blocks. 
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Crossed contribution 



Type 1 



Recall that this contribution describes interference be- 
tween counter-propagating amplitudes. Its graphical rep- 
resentation shown in Fig. [10] is constructed analogously 
to the ladder terms by expanding, on the right hand side, 
each of the single-atom building blocks (C), (F), and (D) 
into a sum of elastic and inelastic responses. 

The diagrammatic expansion of the block (F), 
apparently, contains the new response function 
P' ^(u)\,u)%; v) shown in Fig. Illf a). However, in 
full analogy with the functions (lt ± (cji, lo2))^ \ the 
function P^ ^ (cJi , U2 ; v) is the same as the function 
P( ^(w,u);v) shown in Fig. Hfli), but taken with 



FIG. 10: (Left) single-atom building blocks for the crossed 
contribution of type 1 arising due to interference between the 
counter-propagating amplitudes, see Fig. [5jc). The detected 
signal originates from the blocks (C) and (D); (right) expan- 
sion of the right hand side into the elementary single-atom 
building blocks. 
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FIG. 11: Elementary building blocks required to describe the 
inelastic response of the intermediate atom: (a) block (f5) in 
Fig. [El (b) block (g9) in Fig. \V2\ 
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FIG. 12: (Left) single-atom building blocks for the crossed 
contribution of type 2 arising due to interference between the 
counter- propagating amplitudes, see Fig. E^d). The detected 
signal originates from the blocks (G) and (D); (right) expan- 
sion of the right hand side into the elementary single-atom 
building blocks. 



different frequency arguments corresponding to the 
incoming arrows. 

Combining all contributions (C), (F), and (D) yields 
45 terms. However, 13 of them contain closed loops (see 
Appendix |Pl). and are forbidden according to the rules of 
combining diagrams (see Sec. HVDj) . Only the remaining 
32 triple scattering diagrams (multiplied by a factor 6) 

(3) 

thus contribute to the crossed intensity C tot 1 of type 1, 



see Eq. {22 



2. Type 2 

Expanding the single-atom blocks (A), (G), and (D) 
into the elastic and inelastic responses, we obtain 2, 9, 
and 3 diagrams, respectively, in the right hand side of 
Fig. [T2] Analyzing the diagrams contributing to block 
(G), we notice a new term displayed in Fig. ITTT b). 

Combining single-atom responses on the right hand 
side of Fig. [12] yields 56 terms. Excluding the diagrams 
featuring closed loops (8 diagrams, see Appendix [D]) we 
come up with 46 terms. Finally, the crossed intensity 

(3) 

C tot 2 °f type 2 results as 12 times the real part, see 
Eq. (j2"3")l . of the sum of all these 46 terms. 



FIG. 13: Example of a triple scattering interference spectral 
response obtained by self-consistently combining diagrams 
(a2), (g9), and (dl) from Fig.[H 



D. Self-consistent combination of single-atom 
building blocks 

Figures [51 QJJ] and [TH present diagrammatic ex- 
pansions of the single-atom responses in the elementary 
single-atom building blocks. To obtain the triple scatter- 
ing ladder and crossed spectra from these expansions, we 
combine the elementary blocks in a self-consistent way 

First, we compose the three-atom diagrams using 
single-atom blocks by connecting the outgoing arrows 
with incoming ones regarding the direction and character 
(solid or dashed) thereof. Second, according to the spec- 
tral responses functions associated with the blocks given 
in Figs. [51 El andfTTl we ascribe frequency values to each 
of the inelastic (that is, distinct from the laser frequency) 
arrows in the resulting diagrams. In doing this, we also 
keep in mind that the two outgoing (downward) arrows 
correspond to the spectral intensity of the backscattered 
light at a given frequency, and must have the same fre- 
quency value v. Finally, all intermediate inelastic fre- 
quencies which alter upon scattering are integrated over. 

We already mentioned that 13 diagrams resulting from 
Fig. [10] and 8 - from Fig. [12] are forbidden because they 
contain closed loops wherein the scattered amplitudes 
travel between the atoms without going out (see an ex- 
ample of a loop in Appendix [D]) . We note that similar 
forbidden diagrams appear also in a nonlinear transport 
theory of classical scatters [26( . 

It is instructive to illustrate the above rules by an ex- 
ample. For that, consider a self-consistent combination 
of diagrams (a2)(g9)(dl) (see Fig. IT2"j). which constitutes 
1 of the 46 diagrams contributing to the crossed inten- 
sity of type 2. The corresponding compound diagram is 
drawn in Fig. 1131 The symbolic expressions for the re- 
sponse functions corresponding to the elementary build- 
ing blocks appearing in Fig. [JJJ] are given in Figs. El[a), 
(b), (e), and [TTTb) (more precisely, the block in Fig.^a) 
is the complex conjugate of the one from Fig. [13]) . Not- 
ing that the intermediate inelastic frequency uj changes 
its value, we obtain the following result for the contribu- 
tion: 

/CO 
-OO 

x (a + )(°\a-(v)}(-\ (24) 

where we have omitted the prefactor (T^T^T^T^) conf 
which scales as (fci€) -4 (see Eq. @), where I is the mean 
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.CO, 
V 



7 



<3((D l ,(D 2 ,...,(D n )> 



(+-...-) 



FIG. 14: Circle with incoming dashed and solid arrows at fre- 
quencies u>i, (jJ2, ■ ■ ■ ! Wni representing the nth order correction 
to the optical Bloch vector (<?(u;i, tV2, ■ • • ,0J n )}^ " > whose 
frequencies and signs in the superscript are arranged in the 
same order. The signs ' + ' and ' — ' are associated with dashed 
and solid arrows, respectively. The positive and negative fre- 
quency amplitudes associated with this block have frequency 
SiUJi which is in general non-zero, i.e. detuned from the 
laser frequency, and hence the circle is hatched. 



interatomic distance. Any other combination of the el- 
ementary building blocks is expressed in an analogous 
way. Further examples are considered in Sec. EH 



V. GENERAL EXPRESSIONS FOR 
SINGLE-ATOM BUILDING BLOCKS 

In the previous Sees. IIIII and IIVI we introduced dia- 
grams for all elementary single atom building blocks, see 
Figs. [Jj [U [10] and [J31 exhibiting at most four incom- 
ing arrows. In this section, we will give the correspond- 
ing mathematical expressions, generalized to an arbitrary 
number n of incoming probe fields. 

It is important to have formulae for an arbitrary n for 
the following reasons: First, the laser-driven cold atoms 
in a cloud experience random polychromatic fields which 
can be modeled by a large number of the incident ar- 
rows labeled by different frequencies 



181 ] . Second, al- 



ready in case of three incoming arrows the expressions 
for the building blocks become rather heavy. It is these 
two considerations that urged us to look for a general and 
concise description of the building blocks. Such a formu- 
lation is possible due to the clear hierarchical structure of 
the building blocks and the small number of constituents 
needed for their construction. 



To compose an arbitrary spectral response function, 
we need (i) the single-atom Green's function G(z) — 
(z-M)- 1 , where M = M x according to Eq. (|A2|) with 
Q\ = Q for an atom placed at the coordinate origin; (ii) 
the matrices A^ + ) and A^ - ) describing the coupling of 
the negative- and positive-frequency probe fields to the 
atomic dipole (Eq. (|A8I) ): (iii) frequency values and types 
(dashed or solid) of all the incoming probe fields. 



Spectral response functions associated with the 
Bloch vector 



In the hierarchy of the single-atom building blocks, the 
most basic element is the steady state optical Bloch vec- 
tor for an atom interacting with a laser field alone. We 
will represent this vector graphically by a blank circle: 



O = <<?) (0) = GL, 



(25) 



where L is given in Eq. (TA"2]) . and G = G(0). 

Circles with one, two, n incoming arrows describe 
perturbative corrections to the Bloch vector due to one, 
two,..., n weak classical probe fields. Equivalently, these 
corrections can be regarded as arising from the excita- 
tion transfer processes from the surrounding atoms (see 
Sec. IIII C ip . A graphical representation of the nth or- 
der correction (ct(wi, u>2, ■ ■ ■ ,u n ))^ for a particular 
choice of the incoming arrows is shown in Fig. Q31 

As follows from Fig. [JJJ we associate signs Sj (j = 
l,...,n) with each of the n incoming arrows according 
to the rule 



— , for incoming 
+, for incoming 



(26) 



We stress that, to avoid confusion, we always arrange fre- 
quencies and the corresponding signs in the same order. 

With these preliminary remarks, we intro- 
duce the general expression for the vector 
(a( Wl ,w 2 ,...,w„)}( s i S2 - s "): 



H Wl , W2 ,...,^)>(^- s "> - J2 G (*ELi^J^ J ...G^ (27) 

t(J1i— t3n) 



where . . . , j n ) denotes n! permutations of indices 

ji, . . . , j n € {1, . . . , n}. As can be proven by the method 
of induction, the same expression is obtained when ex- 
panding the quasi-stationary solution of the optical Bloch 
equations for an atom driven by a superposition of the 
laser field plus n additional fields with different frequen- 
cies Ui,...,w n up to first order in each of the additional 



field's amplitudes. 

To illustrate application of the formula ([2"T[) by an ex- 
ample, we will consider the function ((?(wi, W2)/ H ' de- 
scribing a second-order correction to the Bloch vector due 
to the incoming negative and positive frequency ampli- 
tudes at frequencies oj\ and 0J2, respectively (Fig. L^b)). 
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By virtue of Eq. ([27]) we obtain: 

(<t(wi,w 2 ))( + -) = G(i[wi - w 2 ])A(- ) G(iwi)A(+)G J L 
+ G(i[wi - w 2 ])A( + 'G(-ia) 2 )A(-'GL. (28) 

The first and second element of the vector 
..,«„)>(**•"*»), (a^( Wl ,..., W „))^- s ") and 
((T + (o;i, . . . ,( 1 )„))' Sl,,,s °', correspond to the outgoing 
positive and negative frequency amplitude depicted by 
a solid and dashed arrow, respectively. The simplest 
outgoing amplitudes correspond to absent incoming 
arrows: 

O— > = <O (0) =[G2] 1 , (29) 
O- = (<r + )W = [GL] r (30) 

The frequencies of the outgoing arrows are uniquely de- 
fined by those of the incoming ones and are equal to 
— Sj^j for a solid and + • SjUj for a dashed arrow. 
Unless, for each arrow with non-zero detuning ujj there 
exists a corresponding complex conjugate arrow with the 
same detuning, J2j s j UJ j 0- I n Fig. HH we have as- 
sumed that the latter is true, and therefore hatched the 
circle. 



B. Spectral response functions associated with 
correlation functions 

Previously, we described the amplitude response func- 
tions which are depicted graphically by a circle with 
n > incoming and a single outgoing arrows. Now, 
we will consider the building blocks with two outgoing 
arrows, of which one is solid and one dashed. In gen- 
eral, these building blocks are associated with the atomic 



dipole correlation functions 20]. As already mentioned 
above, splitting these functions into a factorizable and 
non-factorizable component then defines the elastic and 
inelastic component of the radiated intensity. 



1. Factorized response functions 

In the beginning, we will consider the building blocks 
which can be factorized as products of the negative and 
positive frequency scattering amplitudes. In the simplest 
case of no incoming arrows, such a product gives the 
intensity of the elastic component of resonance fluores- 
cence. Graphically, this quantity is represented as 

o— > 

(31) 

o— >• 

For an arbitrary number n of the incoming arrows, they 
can be distributed among the outgoing positive and nega- 
tive frequency amplitudes of a factorized correlation func- 
tion in 2™ ways. The resulting correlation function is 
then expanded into a sum of these 2 ra combinations. For 
n < 3, the diagrammatic expansions are given in Figs. [7j 

ma andm 

In the following, we will use permutations between 
the sets of indices {ji, . . . , j k } and {jk+i, ■ ■ -,jn}, which 
W6 will denote by tt(ji , ■ • ■ , jk \ jk-\-i •>■•••> jn 

). There are 

n\/k\(n — k)l such permutations. Now, if there are 
n incoming arrows, the factorized response function 
g( s i— s »)(o; 1) . .. , oj n ) reads 



+ . . . 

+ Yl (^k)) (s,1, (^K.- • • ,Wi»)> ( * ,a -* , - ) 

+ <a+>(°)(a-( Wl ,...,w„))( Sl - s "). (32) 
Let us illustrate Eq. (|32|) by an example for n = 2: 

fl< + ->(wi,U>a) = (<T + {U> 1) W 2 W + - ) (<T-)W + (<7-(wi,W2)> (+ -V+>< > 

+ <tr+(w 1 )>( + )( £ r-(w 2 ))(-5 + (a~ (a + (lj 2 ))^\ (33) 

I 



where the subsequent terms of Eq. (|33]l correspond to diagrams (fl), (f2), (£3), and (f4) with the frequencies of 
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, to, 

\ 2 ■•■ 

( ) 

'• [j^jl . * = P k> 2 ,..., co n ; v) 

' k '** 

FIG. 15: Inelastic building blocks have two outgoing arrows. 
The frequencies of the outgoing and incoming arrows are re- 
lated by the energy conservation relation v — v' = Y.k s k ijJ k 
(see Eq. pi))). 

the incoming dashed and solid arrows equal to oji and 
U)2, respectively (see Fig. ITU1). 

2. Non-factorized response functions 

The fluctuating part of the atomic dipole correlation 
function cannot be factorized. The building blocks as- 
sociated with such functions are depicted by hatched 
boxes with one dashed and one solid outgoing arrows (see 
Fig. [H]), and are denoted by p( Sl — Sn '{u)i, . . . , uj n ; v). In 
this case, the frequencies of the outgoing arrows are not 
defined only by those of the incoming ones: an addi- 
tional frequency v appears which describes the spectral 
distribution of the outgoing, say dashed, arrow. The fre- 
quency v' of the outgoing solid arrow is related to v by 



In case of no incoming arrows, Eq. (|3"6")l reduces to the 
following expressions: 

qf = -iA<+>(*><°> +£i - (a) (0) (O (0) > 



the identity 

71 

v' = v - ^ s k QJk, (34) 

k=l 

which reflects the energy conservation upon the scatter- 
ing process. 

It is well known how to obtain the inelastic part 
of the emission spectrum of a laser-driven atom 
(see, for instance, [HI). To find the function 
p( sl '" Sn ' (u)i, . . . ,uj n ;i/), we have generalized the same 
procedure to account, in addition to the laser field, for 
the presence of n probe fields. 

An important ingredient in the definition of the 
function p( s i— s ")(u;i, . . . ,u; n ; v) are the expressions 
for the fluctuating part of the nth order correction 
of atomic dipole correlation functions. These ex- 
pressions can be straightforwardly obtained with the 
aid of the nth-order correction to the Bloch vector 
(o ! (wi,...,w n )( sl " ,% ', Eq. (p7]l . and the factorized re- 
sponse function g( Sl - Sn ^ (lji, . . . ,w„), Eq. (|32|) . 

Using the function g^ Sl — Sn '(u>i, . . . ,w n ), we cre- 
ate two vector functions, g*+"' (u}i, ■ ■ ■ ,ui n ) and 

g_'" B (oji, ■ ■ ■ ,u tt ). The former and latter vectors are 
obtained by making substitutions er + — >• a and o~ — > <j, 
respectively, in each of the terms of Eq. (|32|) . Explicitly, 



(35) 



(36) 



where the L x = (0, 1/2, 0) T , and L 2 = (1/2, 0, 0) T appear 
due to the identity (cr+cr-) (0) = (1 + (^) (0) )/ 2 - 

Vectors q± S1 ' " Sk \u>i, . . . , u>k) (k = 0, 1, . . . , n) enter the 
definition of the inelastic building block shown in Fig. Q~5] 
through the relations: 



I 

^••• s "Vi, ... = <*(wi, • ■ ,4 (,I "* ) (ir T > (0) 

7rCjli".)Jn-l|in) 
+ . . . 

+ J2 (^• 1 )) Cs « ) (^(^ 2) ... 1 ^j) Csfe --- s - ) 

Now, the fluctuating part of the atomic dipole correlation function reads: 

. ,. )U ,) = T i4(±)(?( Wl , . . . , Wn ))^- S ") - ^-'^(wi, . . .,«„). 

I 



I 

P^- s -)( Wl , . . . ,«„;!/) = ^ (Pi* 1 -- ) (w lj . . .,«„;!/) + P^-'-^wi, . . . , w„; i/)) , (37) 
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where 



r>(si...s„) / \ 



^ [g^A^ . . . A (s ^G(iu' + i s J1 a; J1 )A (s " ) G(^')'f| 0) 
V A^» ) . . . A^ ) + is^ Wil )^J' Jl } (uj h ) 



+ J] [G(-ii/)A('i-)...A^)GHi/ + « il w il )^ l) (a;j i ; 



(38a) 



G(-»i/)^™' , - ) (w il ,...,w iB )' 

I 



(38b) 



with v 1 defined in Eq. (|34p . Again, as can be proven by 
the method of induction, Eqs. (|37H38p coincide with the 
expressions obtained from the perturbative solution for 
a single atom driven by a laser field and additional weak 
probe fields. 

Finally, we will illustrate Eqs. ([37]). ([38]) by the sim- 
plest non-trivial example of the expression for a hatched 
box with one incoming dashed arrow at frequency u> (see, 
e.g., Fig. |T07c3)). In this case we have 



Z7T V L 



(0) 



+ 



G(iu- iv)A^G{-iv)q^ ] 
G(iuj - w)q W (w) 



(39) 



C. Derivation from the master equation 

Having defined the elementary single-atom building 
blocks in Sees . I V Al and IVBl and the rules for connecting 
them with each other in Sec. IIV1 we are now able to cal- 
culate all the various components (ladder and crossed of 
type 1 and 2) of the triple scattering intensity in terms of 
single-atom quantities. It remains to be shown that the 
corresponding expressions are identical to those derived 
from the master equation introduced in Sec. IIIII 

To obtain the solution of the master equation in the 
form of self-consistent combination of single-atom build- 
ing blocks, we generalize the method that was previously 
applied to two atoms [ljj]. Some steps towards such a 
solution have been already made in Sec. IIIII Namely, 
in Sec. 1111 Bl we introduced the recurrence relations for 
vectors x^ n \ y^ n \ and Thereafter, in Sec. IIII Cl 



we selected the triple scattering diagrams for which the 
incoherent and interference intensities survive the disor- 
der average, resulting in the ladder and crossed contribu- 
tions. 

The next step is to factorize all terms containing 
the relevant vectors x ^ , y ^ as combinations of three 
single-atom expressions. Recall that these vectors in- 
clude four interaction matrices sandwiched between the 
Green's functions G + , G x , or D x (see Eq. ([T2"|l ). The 
factorization is then accomplished by using the integral 
representations for G x and D x which express the two- 
and three-atom evolution matrices as frequency integrals 
over tensor products of the single-atom Green's matrices 
(see Appendix iBj) . 

After this procedure, each term contributing to the lad- 
der and crossed intensity contains multiple (up to seven- 
fold) frequency integrals. Although each integrand repre- 
sents a product of three single-atom expressions, neither 
of them can be interpreted physically. To obtain the 
physically meaningful single-atom building blocks as de- 
fined by the solution of the optical Bloch equations for a 
single atom, we simplify these expressions further. 

To this end, by using the sum rules (see Appendix [Cj 
we solve exactly some of the frequency integrals. As a 
result, we obtain the final expressions for the ladder and 
crossed spectra which include at most two integrations 
(see, e.g., Eq. (|4^|) valid for very strong driving £1 3> 7). 
The ladder spectrum (Fig. [5ja) and (b)) is given by a 
sum of 118 = 50 + 68 terms corresponding to all possible 
combinations of the single-atom building blocks shown 
in Fig. [ 3(50 = 2 x 5 x 5) and Fig. M (68 = 2 x 17 x 2), 
see Sec. IIVB II above. Likewise, the crossed spectrum 
(Fig.[5jc),(d)) is given by all allowed combinations of the 
single-atom building blocks shown in Figs. [TU] and fSl 
see Sec. IIV 01 This completes the derivation of the triple 
scattering contribution expressed in terms of single-atom 
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FIG. 16: Triple scattering diagrams yielding the elastic (a-e) 
and inelastic (f-h) ladder spectra of type 1. For one-photon 
processes (~ fi 2 ), only diagram (a) contributes. 



response functions. 



VI. TRIPLE SCATTERING SPECTRA 

In this section, we present our results for the triple 
scattering spectra of CBS. In the beginning, we provide 
analytical expressions for the elastic and inelastic con- 
tributions in the limit of small Rabi frequency <C 7. 
Thereafter, we consider the non-perturbative regime of 
strong atom-laser interaction. 



A. Case of small Rabi frequencies (Q, -C 7) 

In this limit, it is natural to restrict ourselves to the 
triple scattering contributions that are not larger than ~ 
(f2/7) 4 . Physically, this corresponds to a weakly inelastic 
scattering by three atoms of the laser field containing not 
more than two photons. 

Analyzing the behavior of the triple scattering dia- 
grams, we established that they exhibit different asymp- 
totics when ft -C 7. For example, the crossed diagram 
(a2)(g9)(dl) (see Eq. ([211)) yields an expression which is 
~ (fi/7) 6 , and can therefore be neglected. We will con- 
sider separately the elastic and inelastic components of 
CBS. 



Elastic spectrum 



Diagrams in Fig. 
tic ladder intensity 



16Ta)-(e) describe 



-(3) 



They 



the elas- 
be pre- 



sented in the form of equations: (a) = (al)(b3)(b3), 

(b)=(al)(b3)(b2), (c)=(al)(b3)(bl), (d)=(al)(bl)(b3), 

and (e)=(al)(b2)(b3), where the compounds on the 

right hand sides are the single-atom building blocks 

from Fig. Figure I17f a-d) depicts diagrams contribut- 
or 

ing to L ei2 - In this case, the corresponding equa- 
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FIG. 17: Triple scattering diagrams yielding the elastic (a-d) 
and inelastic (e) ladder spectra of type 2. 
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FIG. 18: Triple scattering diagrams yielding the elastic (a-e) 
and inelastic (f-h) crossed spectra of type 1. For one-photon 
processes, only diagram (a) which is a reciprocal of diagram 
(a) in Fig. 1161 contributes. 



tions for diagrams read (see Fig. [SJ: (a)=(al)(e7)(al), 

(b) =(al)(el3)(al), (c)=(al)(e6)(al), (d)=(al)(el6)(al). 
Likewise, the elastic crossed contributions in 

Fig. [T8l a)-fe) result from the combinations of single-atom 
blocks of Fig. [TJ3 (a) = (cl)(f3)(dl), (b)=(cl)(fl)(dl), 

(c) =(cl)(£2)(dl), (d)=(c2)(f3)(dl), (e)=(cl)(f3)(d2), 

(3) 

which yield the elastic crossed intensity C cl \ . Finally, 
the four diagrams in Fig. [19] yield the crossed contri- 
bution Cgfg- They are composed of the diagrams of 
Fig. [H as follows: (a) = (al)(gl)(dl), (b)=(al)(g2)(dl), 
(c)=(al)(g5)(dl), (d)=(al)(g7)(dl). 

We obtain the following analytical expressions for the 
elastic triple scattering ladder and crossed intensities 
which are valid up to (f2/7) 4 : 



r(3) 



2m 4 



r(3) _ 

— 



r (3) _ 
°el,l — 

(3) _ 



32(1 + <5 2 ) 3 32(1 + 5 2 ) 4 ' 
6ft 4 



32(1 + 6 2 ) 4 ' 

3ft 2 



24ft 4 



32(1 + 6 2 ) 3 32(1 + <5 2 ) 4 ' 
24ft 4 



cl,2 



32(1 + S 2 ) 4 



(40a) 
(40b) 
(40c) 
(40d) 



where ft — ft/7, 5 = S/j. Using the above results, 
which, as we have checked, coincide with the results 
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FIG. 19: Triple scattering diagrams yielding the elastic (a-d) 
and inelastic (e) crossed spectra of type 2. 



following from the two-photon scattering theory we 
draw two conclusions. First, in the linear regime, only 
the diagrams in Figs. fTBTa) and ITST a) contribute to the 
signal, and ladder and crossed components are equal to 
each other, what indicates perfect phase coherence in the 
linear regime. From Eqs. (|40a[) and (|40cl) it follows that 



r(3) _ ^(3) 
-^etl — °cl,l 



3ft 2 



32(1 + S 2 ) 3 



(41) 



Second, at increased Rabi frequencies, the interference 

(3) (3) 

contribution , + C cl 2 decreases faster than the ladder 



contribution L cll + L^ J 2 , since the contributions propor- 
tional to Cl 4 have a negative sign, and are larger for the 
crossed component (24 + 24 = 48) than for the ladder 
component (27 + 6 = 33), see Eqs. ([i0|. The fact that 
the interference contribution exceeds the background can 
be explained by the fact that, in the case of nonlinear 
triple scattering, the coherent backscattering signal is in 
general formed by three interfering amplitudes 9] . 



(3) 



2. Inelastic spectrum 

The four inelastic ladder diagrams are depicted in 
Fig. ll6f fV fh) and Fig. [T7T e). These diagrams result from 
the combinations: (f)=(al)(b3)(b5), (g) = (al)(b5)(b3), 
(h)=(a2)(b3)(b3) (see Fig. EJ), (e) = (al)(el7)(al) (see 
Fig.©. 

Concerning the inelastic crossed spectrum, it is given 
by the four contributions shown in Fig. [TBT fl-fh) and 
Fig. [TWcV Diagrams (f)-(h) consist of the following 
building blocks: (f)=(cl)(E)(dl), (g) = (cl)(f3)(d3), and 
(h)=(c3)(f3)(dl), with the right hand sides coming from 
Fig. [TO] "Type 2" crossed diagram in Fig. ITW e) yields 
one contribution to the inelastic spectrum in the weakly 
inelastic regime: (e) = (al)(g9)(dl), with the constituents 
defined in Fig. [T^l Due to this contribution which 
emerges when the scatterers are nonlinear 0, [26[ , the in- 
terference effect modifies dramatically: the enhancement 
factor associated with the inelastic scattering becomes 
greater than 2, again due to the interference between 
three amplitudes mentioned above. 



xlO 




FIG. 20: (Color online) Inelastic ladder and crossed spectra 
of triple scattering for the case of exact resonance (8 = 0) and 
Q = O.OI7. 



Indeed, consider the analytical expressions for the in- 
elastic ladder and crossed spectra: 



4 3 L(*) 



1 



m 4 



2tt 32(1 + <5 2 ) 3 
(3(1 + 5 2 ) + 45v - 



(42a) 



1v 2 f 



L 



(3) 



incl.2 



(1 + (5 - i>) 2 )(l + (5 + 9) 2 f 

1 n 4 



2tt 32(1 + S 2 ) 3 
12 



(l + (8-i>) 2 )(l+{S + i>) 2 ) 

1 6n 4 



d 3) 



incl,2 



2tt4(1 + ^ 2 ) 3 

(1 + S(6 + D)) 2 
(l + (8-i>) 2 )(l + (8 + P) 2 ) 3 ' 

1 6n 4 



(42b) 



(42c) 



2tt4(1 + £ 2 ) 3 

1 + 8(8- 



(42d) 



(1 + (5 - u) 2 )(l + (5 + v) 2 ) 2 

As we have checked, the inelastic spectra given by 
Eqs. (|42|) are also in full agreement with the two-photon 
diagrammatic scattering theory Q. The curves corre- 
sponding to the functions L^u) = ij^iM+lj^jM 



CinliM + & a M. are plotted in 



t(3) 



and Cgiiv) 

Fig. m 

It is clearly seen that the dashed line corresponding to 

Qnci(^) 1S above the solid line corresponding to L\J cX (v). 
To characterize this effect in more quantitative terms, let 
us, following Q, assume that the elastic component has 
been filtered out, and study the behavior of the enhance- 
ment factor, which in this case is defined as 



C, 



(3) 
incl 



(3) 



(43) 
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1. Elastic spectrum 




FIG. 21: (Color online) Elastic intensities as a function of 
Rabi frequency at exact resonance: (solid) (dashed) 



C. 



(3) 



Inset displays two elastic crossed contributions: 



(dashed-dotted) C$\; (dotted) C, 



-.(3) 



The elastic intensity exhibits an oscillatory-like behav- 
ior as a function of the Rabi frequency for < 27. The 
intensity can assume negative values, which is in accor- 
dance with the fact that the triple scattering is observed 
from the laser-driven transition. Therefore, it is only a 
correction to the physical intensity which also includes 
single scattering, and is always positive. 

A curious feature exhibited by the crossed contribu- 
tions in Fig. |2"TT inset) is that, in a range of Rabi frequen- 
cies where the elastic intensity is significant, they have 
opposite phases. For values of < 0.5j, C ol 2 < < 

CV r If O.57 < Q, < 7, both crossed contributions fea- 
ture destructive interference. In the range 7 < Q < 2j 
the crossed contributions again exhibit opposite interfer- 

(3) (3) 

cncc character: C£{ < < C£ J 2 . 

Although we have not found a simple explanation of 
why the two crossed contributions exhibit opposite inter- 
ference character, this fact is not surprising given that 
they originate from different triple scattering paths. 



Integrating Eqs. (|4"2"]l over D, we obtain: 
, (3 ) _ 3(22 + ~6 2 )& 



/-TV 



L 



128(1 + 5 2 ) A 
(3) _ 3(154 + 25£ 2 + 3£ 4 )» 4 
~ 1024(1 + 1 2 ) 4 : 



(44) 
(45) 



wherefrom we deduce that r](6) > 2 for \6\ < 1.042, with 
the maximum value of ~ 2.143 at 8 = 0. In fact, in a 
cloud consisting of a large number of atoms, and applying 
appropriate frequency filtering, the enhancement factor 
(ffS)) can in principle reach the value 3 [23| . 

So far, we have restricted our consideration to the case 
of weakly inelastic scattering which is valid for very small 
Rabi frequencies. We will next address the behavior of 
the CBS spectra in the non-perturbative regime of the 
atom-light interaction. 



B. Non-perturbative account of atom-laser 
interaction 



If the Rabi frequency does not satisfy the condition 
Q -C 7, we enter the regime of light scattering by atoms 
wherein the non-perturbative expressions for the single- 
atom building blocks must be used. Consequently, we 
consider contributions from all triple scattering diagrams 
resulting from self-consistent combinations of the single- 
atom building blocks shown in Figs. [3 13 [101 an d HU 
The results of our numerical calculations of the triple 
scattering elastic and inelastic spectra are presented in 
Figs. [5l] and [2U respectively. 



2. Inelastic spectrum 

To see how inelastic processes affect CBS from three 
atoms, in Fig. [22] we plot the ladder and crossed spectra 
for different values of the Rabi frequency. For ft = O.I7, 
the two-photon processes still give the dominant contri- 
bution to triple scattering, which is manifest in a qual- 
itative agreement between the spectra in Figs l2"27 a) and 
[20l However, at increased values of the Rabi frequency, 
Fig- l2"27 b-f). the crossed spectra are dominated by the 
ladder spectra. Furthermore, small in magnutude inter- 
ference spectra feature destructive interference, or CBS 
anti-enhancement (see Fig. [227 c-f)). due to the crossed 
contribution of type 2 (see insets). 

Increasing the Rabi frequency is also accompanied by 
a splitting of the ladder and crossed spectra into several 
Lorentzian and dispersive resonances. This transforma- 
tion of the spectra can be attributed to the dressing of 
atomic states by the laser field [16J. 

For very strong driving strengths 3> 7, the in- 
tensities of the elastic components become vanishingly 
small, and the total spectral line shapes are given by the 
corresponding inelastic spectra. The latter can be de- 
duced by combining the three elementary inelastic build- 
ing blocks (hatched boxes) for each of the contributions. 
Namely, the combinations (a2)(b5)(b5), (a2)(el7)(a2), 

(c3)(f5)(d3), and (a2)(g9)(d3) yield A 2 {v) t 



and C\ 



C 



and [T 



inii.aM' respectively (see Figs. H HOI 
The inelastic ladder and crossed spectra can 
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FIG. 22: (Color online) Inelastic ladder (solid) and crossed 
(dashed) spectra of triple scattering for the case of exact res- 
onance (5 — 0) for different Rabi frequencies: Q = (a) O.I7; 
(b) O.57; (c) 7; (d) 37; (e) 57; (f) IO7. Insets show two 
inelastic crossed contributions Cf^i 1 (p) (dashed-dotted) and 
c£i 2 (*) (dotted). 



then be presented by the compact formulas: 

L £L,i(y) = 6 / dwW'p(°V) (46a) 

x p(+~) (a/, J- W ")P(+-'( W ", oj"; v), 
L£l 12 (p) = 3 / cL/^"P (0 V)-P ( °V') (46b) 



C^^i/) =6 / dw'dw"pW(w / ;i/)P^) 



(46c) 



CtoeiaM = 12Re / dwWP (0 V)-P ( ~V» 



tions and spectral decompositions of the Green's matri- 
ces, as well as the sum rules were developed and applied 
to prove the equivalence between the master equation and 
the pump-probe approaches for double scattering. Here, 
we have generalized these techniques for three atoms. 

From the rigorous analytical solution of the master 
equation in quadratures, we have found the explicit 
analytical formulas for the single-atom building blocks 
whose diagrammatic self-consistent combination yields 
the triple scattering ladder and crossed spectral intensi- 
ties. The obtained expressions for the single-atom spec- 
tral response functions are equivalent to the pump-probe 
solutions following from the optical Bloch equations un- 
der classical bi- and trichromatic driving fields. 

In contrast to double scattering, triple scattering fea- 
tures two crossed contributions in the inelastic scattering 
regime which manifest the multi-wave character of CBS 
[26j |. It is the interplay between these two contributions 
that defines the behavior of the overall interference ef- 
fect as a function of the Rabi frequency. In the weakly 
inelastic scattering regime, the two inelastic crossed con- 
tributions are positive and in phase, yielding the CBS 
enhancement factor > 2 - in qualitative agreement with 
the earlier prediction for nonlinear classical scatterers Q . 
However, increasing the laser intensity leads to a rapid 
decrease of the relative intensities of the crossed contri- 
butions with respect to the ladder term, on the one hand, 
and to the destructive interference, on the other hand. 

The results of this work go beyond the problem of 
triple scattering. We have deduced general analytical ex- 
pressions for single-atom building blocks for an arbitrary 
number of incoming probe fields. These expressions can 
be incorporated into the stochastic pump-probe approach 
[l7L [lH for an accurate simulation of the multiple inelas- 
tic scattering of laser light in cold atomic ensembles. 

Another direction of future work will be to generalize 
our present results for realistic atoms probed in the CBS 
"saturation" experiments This can be achieved 

by including the polarization degree of freedom into the 
description of the single-atom response functions. 



x P(+-+V,u;>-u/». 



(46d) 



VII. DISCUSSION AND CONCLUSION 



We have studied the coherent backscattering of intense 
laser light from three two-level atoms using the master 
equation approach. The goal of this work was to find the 
analytical solution of the problem and show its equiva- 
lence with the solution following from the diagrammatic 
pump-probe approach to CBS. Wc have fulfilled this task 
in the spirit of the earlier work [19[ , where technical tools 
such as the recurrence relations, the integral representa- 
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Appendix A: Evolution matrices 
1. Matrix A: evolution of independent atoms 

The blocks of matrix A (see Eq. (fTTjl ) read 

M + = M 3 © M 2 © Ml, (Ala) 
M x = (M 2 © 1 + 1 © M 3 ) © (Ml © 1 + 1 ® M 3 ) 
© (Mi © 1 + 1 © M 2 ), (Alb) 

(L® 1 1 © L \ 
L© 1 1®L , (Ale) 
L® 1 

L x = (£®1®1 1©L©1 1®1®L), 

(Aid) 

7V X = 1 © 1 ® M 3 + 1 © Mi ® 1 + Mi © 1 © 1. 

(Ale) 

Here, 1 denotes the unit 3x3 matrix, M\ (A = 1, 2, 3) 
is a position-dependent 3x3 matrix describing the evo- 
lution of the Bloch vector of atom A, and L is a three 
dimensional vector 19]. Explicitly, 

/-j + iS o -%a x /2\ 

M A = - 7 - it iQl/2 1 , L = (0,0, -2 7 ) T 

\ -inj -2 7 J 

(A2) 

To be self-contained, we remind here the form of the op- 
tical Bloch equation: 



(a x ) = M x (a x ) + L. 



(A3) 



Given the vector L, we can also define the vector A (see 
Eq. ®): 



A = (L,L,L,0) T , 

where the null vector in (|A4[) contains 54 zeroes. 
Using Eqs. (|Ala[) and (Tl4l) we obtain 



(A4) 



G + (z) = G 3 (z)®G 2 (z)®G 1 {z), (A5) 
with G A (z) = (z - Mx)- 1 , cf. Eq. CEU 



2. Matrix 1/: evolution of the dipole-dipole 
interacting atoms 

A convenient way to specify the blocks {U Xll ) a , (WajiJo 
(a =-i,i_, x) of the matrix V, see Eqs. (fTT|) and ([15]) . 
is implicit - by describing their action on test vectors. 
This is how it was done in the case of two atoms [lj|. 
In fact, due to the pairwise nature of the dipole-dipole 
interaction, the results for two atoms can, after some 
modifications, be adapted to the iV-atom case. Here, we 
will sketch this procedure for three atoms. 



To this end, we will for the moment abstract ourselves 
from the presence of a third atom, and consider how the 
blocks of the 'contracted' dipole-dipole interaction matrix 
(denoted as v) act in the two-atom case. When this is 
done, it is straightforward to define the action of the 
blocks of the matrix V. 

Let ax and a M be 3-component column vectors, with 
the indices A, /i referring to the atom's number. Assum- 
ing A < (j,, we create two test vectors 

= (oju, ax) T , {pxp)2 = a x © a M . (A6) 

The 'contracted' interaction matrices v a (a =-i, i_, x) 
are characterized by the following identities [19| : 

(v^HSx © *„) = ( ^ (A?&) 

(«^)-.(3a ® *„) = ( 2iT ^ A (%z xKh ) . ( A7b ) 



K m )^(«a ©a^) 



2iT A * M A(-)a A [a^]i ^ ' 



(vA M )x(«A®a M ) = -2T Aai A( >a A © A( + 'a t 



(A7c) 



(A7d) 



(A7e) 



(w M A)x(aA©a M ) = -2T MA A (+) a A ©A ( - ) a M , (A7f) 
(«v)x(oa © Up) = -2T^A^ax © A(+»a p , (A7g) 
(^A)x(aA©a M ) = -2T; A A(+)a A © A^a M , 



VA; 



«A 

a,, 

«A 

ax 





(A7h) 


rli © (2 i T A , l A (+) a, I ), 


(A7i) 


2zT MA A(+)a A ©ni, 


(A7j) 


(-2iT A * M A^)a A )©n 2 , 


(A7k) 


n 2 © (-2iT; A A<-)a Al ), 


(A71) 



where [3x]k and [a^k refer to the k-th component 
(k = 1,2,3) in the basis defined by the single-atom 
Bloch vector (a x ) = {(a^), <a+), (af)) T , n x = (±,0,0) T , 
n 2 = (0, i, 0) T , = (0, 0, 0) T , and 








-i/2\ 
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Identities (|A7[> show that the blocks of the matrix w ac- 
complish mutual mappings between the vectors of the 
form pi and p 2 (we omit for the moment the indices A, 

vi ■ Pi -> Pi', «x : p 2 -> P 2 , w L : pi -» p 2 \ (A9) 
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where p- (j — 1,2) is the result of the transformation 
which has the same structure as the vector pj . 

Returning to the three-atom case, it can be shown that 
the matrices U a and W a {a =-i, l, x), acting on the test 



vectors 

Pi = {a 3 ,a 2 ,a 1 ) T , (AlOa) 

P 2 = {3,2 ® S3, oi <g> 03, oi <g> a 2 ) T , (AlOb) 

P 3 = ai ®a 2 ® a 3 , (AlOc) 



map them to each other according to the transformation 
rules: 

C/n :P 2 -> P/, C/ L : Pi -> P 2 ', U X :P 2 ^ P 2 ' , 

(Alia) 

:P 3 ->■ P 2 ', W L : Pi -> P 3 ', VK X : P 3 -> P 3 '. 

(Allb) 

Now, the action of the matrices {U\^) a , (W\n) a , 
(A, n = 1,2,3) on the vectors Pj is fully determined by 
the action of the matrices {v\^,) a on the vectors {p\^)i 
and {px f _ l ) 2 . In other words, the explicit form of the vec- 
tors P-' can obtained by appropriately constructing the 
three- atom vectors out of the vectors {p\^)i and {p\S)2, 
the latter vectors being defined in the right hand side of 
Eqs. (TA71) . 

Vector P[ is created from the vector {flu) 1 ky append- 
ing the null vector at the position corresponding to 
atom (3 =/= X ^= fi, such that the atomic indices of the re- 
sulting vector are arranged in decreasing order, as in the 
vector Pi. Vector P 2 is created from the vector {p x ' fl ) 2 
by appending the 2 vectors ® 0, such that the atomic 
indices are arranged in the same order as in vector P 2 . 
Finally, the vector P 3 is created form the vector {p{^) 2 by 
appending a tensor product with vector ap (/3 7^ A 7^ /x), 
such that the atomic indices are arranged in the same 
order as in the vector P 3 . 

Appendix B: Integral representations for the 
Green's functions 

From the definitions of Eqs. (|Alb[) and (|T1]). it follows: 

G x {z) = G 23 {z) ® G 13 {z) ® G 12 {z). (Bl) 

Here, Gx^z) is the Green's matrix for two independent 
atoms A and fj,: 

G^{z) = {z - M^y\ (B2) 

where is the evolution matrix of two free atoms, 

M A/1 = M A ® 1 + 1 <8> M u . (B3) 

As proven in Ref. 19], the Green's matrix G\^{z) can 
be represented as 

Gx,{z) = J~ ^Ga (I ± «/) ®G M (f T to;') , (B4) 



where Ga(z) is defined in Eq. (fT4|) . Hence, G x (z) (cf. 
Eq. (JbTJ) ) is given by a direct sum of three integrals in 
the form of Eq. ([B4| . 

Using this result, we will show that a similar integral 
representation can be obtained for the Green's matrix 

D x {z) = {z-N x )- 1 , (B5) 

governing evolution of three independent atoms. Matrix 
N x , specified in Eq. (|Ale[) , can be written as 

N x =Mi®1®1 + 1® M 23 . (B6) 

Obviously, N x has the same structure as the matrix M x , 
see Eq. (|B3j) . Therefore, we can represent it, by analogy 
with Eq. dBH, as 

D x (z) = |°° (| ± ij) ®G 23 (| T to/) . (B7) 

Finally, using the integral representation for the Green's 
function G 23 (z'): 




(B8) 

where z' = z/2 =p iu/, we obtain 

(SGaf^iiw'^gGafyTw"^, (B9) 

which yields the integral representation for the function 
D x {z). 

Appendix C: Sum rules 

The following sum rules constitute a method of calcu- 
lating exactly sums of the frequency integrals appearing 
in the expressions for the vectors x^ n \ y ^ , and z ^ 
due to the integral representations of the matrices G x {z) 
and D x {z) given by Eqs. (|B4[) and (|F39|) . respectively. 

The idea of the method is based on the spectral de- 
composition of the single-atom Green's matrices G\{z), 
followed by calculation of the frequency integrals in the 
complex plane and, finally, a summation of the results. 
We will outline here the main steps in deriving the sum 
rules, and illustrate them with two characteristic exam- 
ples. Further useful, though restricted to the two-atom 
case, examples can be found in [19] . 

We will first note that the Bloch matrix Ma govern- 
ing the evolution of atom A (see Eq. (|A2[1 ) allows for a 
spectral decomposition 

M A = ]TrMi A , (CI) 

k 
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where 



l"fc)("fcl 
<«kl«fc> 



A,- 



(C2) 



is the projector on the subspace corresponding to the 
eigenvalue rfc, and \u%) {\v%)) are the right (left) eigenvec- 
tors of the matrix M\ satisfying (vfc \u%,) = for k ^ k' 
[27| . The eigenvalues of the matrix M\ are the roots 
of the characteristic polynomial thereof: 

f(z) = (z + 2 7 )((z + 7 ) 2 + 6 2 ) + (z + 7 )|tt| 2 . (C3) 

While the explicit form of rfc will be unimportant to us, 
we will rely on the fact that for arbitrary |f2| and 5, the 
roots are non-degenerate, and Re [rfc] < 0. 

To illustrate the sum rules, we will calculate the zeroth- 
order correlation functions y W and z (°> . Although the 
results are known without any calculations, cf. Eq. (|13p . 
we will obtain them alternatively by using the sum rules. 
These derivations highlight all steps needed for tackling 
more complicated cases. 

Taking n = in Eq. (fl2b)l . we obtain 



V 



(C4) 



which follows directly from Eqs. (jCll) and (|T4J) , we arrive 
at 



duj' 



2tt 



^—^ \ iu> — rfc 



1 



-iui' - n -ri ' 



-L 



1 



iw' - r fc -r k 



-Li 



1 



-lUj' — ?*; 



-L 



(C8) 



This integral can now easily be taken using the residues 
theorem. By noting that, in the upper complex half- 
plane of uj' ', the integrand of (|CJ8|) has simple poles at 
lu' = —irk, we obtain the result 



G\n = 



fej V 



1 - P A - P M 
L -L+±!^L®—1 

-r-fc - n -n -r fe -rfc - r ; 



1 



V — ® -^-i = G A L <g> G„L, 
^ -r fe 



a.'.; 



(C9) 



Using Eqs. (jAlcl) . ([BTI) . and (|B4|) . we can bring Eq. JC4 
to the form 



whence, in agreement with Eq. (|13b|) . it follows 



y 



(0) 



{G23, Gl3,Gl 2 ) 



(C5) 



where 



Gx» = % {g x (ilu')L ® G M (-iu')G ll L 



\-Gx(iu')G x L<S)G^-iuj')L 



(C6) 



y (0) - {G 2 L®G 3 L,G 1 L(g>G 3 L,G 1 L®G 2 L) T . (CIO) 



We will proceed with a calculation of z ^ for which 
Eq. (Pc) yields 



and, for definiteness, the upper sign in the integrals (|B4|) 
was taken. 

Next, inserting into (|C6[) a spectral decomposition of 
the single-atom Green's matrix 



G x (z)=J2 P " 



Z-Tk 



(C7) 



(0) 



(Cll) 



Inserting into (|C11|) the integral representation (|B9|) of 
D x = D x (0) the definition of L x (Eq. (|Ald|) ), as well 
as the result (|C10|) . we get 



' v ' dw' duo" 
~2 



fGi(iw')^ ® G 2 (ioj'i)G 2 L <g> G 3 (-iw£ )G 3 i + Gi(iu/)Gi£ ® G 2 (iuj'l)L ® G 8 (-tw+)G 3 .L 

7T 27T V 

G 1 (z W ')GiL®G 2 (z^)G 2 L®G 3 (^)L) = / / ^l^Ls 123 (a/, w "), (G12) 



with = ± a//2. By virtue of the spectral decomposition (|C7I) the vector function Qi 23 (uj' , uj") yields the 
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expression 



/pi i p2 i p3 

&23(w , ,w")= E ■ T, no — L ® 7, ,/o 4^ 

* — ' » zo;' — r k iuj" — iw ''1 2 — r; — r; — zcj" — ioj / 2 — r m — A m 



k.Lm 



1 P 1 _ P 2 1 P 3 - 

+ —^-L (g) P- L ® ^- —^-L 

iu>' — Tk —r k iu" — ilo' jl — n —ius" — ilo' jl 



1 m ' m 



1 P 1 -. 1 P 2 -, P 3 -A 

+ —L <g> ; — L ® 2^ L ] . (C13) 

loj' — r k — r k iw" — ltd' /2 — r\ —r\ —iu>" — iuj'/2 — r m J 

Taking integral over u" in the complex plane, we observe that each of the lines of (|C13j) has a simple pole u" = 
—iri + u//2 in the upper half-plane (and, also, oj" = ir m + u//2 - in the lower half-plane), which gives 

oo 27r Vzo;' - r k -77 -iuj 1 - n - r m -r m 

1 p lr^ p i* 1 p~ 



+ — ^L®-^L' 

iui' - r k -rk 1 



ioj' - r k -rk 1 —iw - ri - r m -r m 

1 P 1 -> P 2 -. P 3 -A 
+ — . ; m L). (C14) 



Finally, the integral over u) 1 yields 

/ pi p2 i p3 pi p2 i p3 

f (o) = y- £i^ii_^ ± S^Z+-^L®^L® -Zz-L 

rr 1 \ 1 ~n ~rk -n — r m -r m -r k 1 -r k - n - r m —r m 

ft,t,m 

+fLt ® ® ^2 L^j = V -^LL <» ® — = GiL ® G 2 L ® G 3 L, (C15) 

-rk —n -r k - n - r m J f-^ -r k -r; -r m 

which confirms Eq. (|13c[) . 



It should be mentioned that the sum rules are valid 
also in a more general case when the detunings and/or 
Rabi frequencies depend upon the atomic index A. In- 
deed, the fact that the eigenvalues r£ in these cases de- 
pend on A changes nothing in the preceding derivations. 
That means that the results obtained in this paper can 
be generalized to include, e.g., residual thermal motion 
of atoms via individual Doppler shifts and, hence, de- 
tunings 6\ , or wave propagation in the effective medium, 
in which case the laser field strength (Rabi frequency) 
would depend on the atomic position inside the medium. 

Appendix D: Forbidden diagrams 

1. Combinations of diagrams from Fig. HOl 

The following combinations are forbidden: 
(cl)(fl)(d2), (cl)(f4)(d2), (c2)(fl)(d2), (c2)(f2)(dl), 



(c2)(f2)(d2), (c2)(f2)(d3), (c2)(f3)(d2), (c2)(f4)(dl), 

(c2)(f4)(d2), (c2)(f4)(d3), (c2)(f5)(d2), (c3)(fl)(d2), 

(c3)(f4)(d2). All these diagrams as well as the ones 

listed below contain closed loops (see an example in 

pig. Eg. 

o^o^o 

xx x 

o<---o o 

v i 

FIG. 23: (Color online) Diagram (cl)(fl)(d2) contains an am- 
plitude that comes in and cycles between the circles (high- 
lighted with red) without the outgoing arrow. 



2. Combinations of diagrams from Fig. 1121 

(al)(g4)(d2), (al)(g5)(d2), (al)(g6)(d2), (al)(g8)(d2), 
(a2)(g4)(d2), (a2)(g5)(d2), (a2)(g6)(d2), (a2)(g8)(d2). 
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